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Abstract 

^vq " This paper is an updated and extended version of our note [TJ (cf. also [2]). 

To compute difference Grobner bases of ideals generated by linear polynomials we 
adopt to difference polynomial rings the involutive algorithm based on Janet-like 
division. The algorithm has been implemented in Maple in the form of the package 
LDA (Linear Difference Algebra) and we describe the main features of the package. 
Its applications are illustrated by generation of finite difference approximations to 
linear partial differential equations and by reduction of Feynman integrals. We also 
present the algorithm for an ideal generated by a finite set of nonlinear difference 
polynomials. If the algorithm terminates, then it constructs a Grobner basis of the 
ideal. 
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1 Introduction 
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Being invented 47 years ago by Buchberger [3] for algorithmic solving of the membership 
problem in the theory of polynomial ideals, the Grobner basis method has become a 
powerful universal algorithmic tool for solving various mathematical problems arising in 
science and engineering. 

Though the overwhelming majority of Grobner basis applications is still found in com- 
mutative polynomial algebra, over the last two decades a substantial progress has also 
been achieved in applications of Grobner bases to noncommutative polynomial algebra, to 
algebra of differential operators and to linear partial differential equations (cf., for exam- 
^ pie, the book 0]). As to the difference algebra, i.e. algebra of difference polynomials [5], 
in spite of its conceptual algorithmic similarity to differential algebra, only a few efforts 
have been made to extend the theory of Grobner bases to difference algebra and to exploit 
their algorithmic power [3 El [7] . 

Recently, three promising applications of difference Grobner bases were revealed: 

• Generation of finite difference approximations to PDEs [H [9] . 

• Consistency analysis of such approximations [TUl [TT]. 

• Reduction of multiloop Feynman integrals to the minimal set of basis integrals [12] . 

In this note we describe an algorithm (Section H]) for constructing Grobner bases for linear 
difference systems that is an adaptation of the polynomial algorithm [T3] to linear differ- 
ence ideals. In so doing, we construct a Grobner basis in its Janet-like form (Section [3]), 
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since this approach has shown its computational efficiency in the polynomial case [131 [H]. 
We briefly outline these efficiency issues in Section [5j The difference form of the algo- 
rithm exploits some basic notions and concepts of difference algebra (Section [2]) as well 
as the definition of Janet-like Grobner bases and Janet-like reductions together with the 
algorithmic characterization of Janet-like bases (Section [3]). Extension of the notion of 
Grobner basis to nonlinear difference polynomials, which has not been addressed in PQ, [2], 
is briefly described in Section |6] where we also present the algorithm [TJJ for construction 
of such bases. In Section [7] we present our Maple package LDA for computing Grobner 
bases of linear difference ideals, i.e. ideals generated by linear difference polynomials. The 
package is a modified version of our earlier package [T7] oriented towards commutative 
and linear differential algebra and based on the involutive basis algorithm [TJ] . The mod- 
ified version is specialized to linear difference ideals and uses both Janet and Janet-like 
divisions [13] adopted to linear difference polynomials [15] . In Sections |8] and [9] we illus- 
trate LDA by simple examples of its application to the construction of finite difference 
approximations to linear systems of PDEs and to the reduction of Feynman integrals. 

2 Elements of difference algebra 

Let {y 1 , . . . , y m } be the set of indeterminates, e.g., m functions of n variables xi, . . . , x n , 
and {6*1, . . . , 9 n } be the set of mutually commuting difference operators (differences), i.e., 

(6i o y J )(x!, ...,x n ) = y J (x ly . . . ,x { + 1, . . . ,x n ). 

A difference ring R with differences 9i, . . . , 9 n is a commutative ring R such that for all 
f,g e R, 1 < i,j < n, 

CO, OjO;. e i o(f + g)=9 i of + e i og, 
i o(fg) = (0 l of)(9 l og). 

Similarly, one defines a difference field. 

Let K be a difference field, and 1Z := K{y 1 , . . . , y m } be the difference ring of polyno- 
mials over K, in variables 

{0»oy k \ fieZ£ Q ,k = l,...,m}. 

Hereafter, we denote by IZl the set of linear polynomials in TZ and use the notations: 

= {0" | ^eZ| }, deg l (^oy fc )= /ii , 
deg(^ o y k ) = \n\ = J2i =1 Hi ■ 

A difference ideal is an ideal X C TZ closed under the action of any operator from 0. 
For F C TZ, the smallest difference ideal containing F will be denoted by Id(F). If for an 
ideal X there is F C TZl such that X = ld(F), then J is a linear difference ideal. 

A total ordering y on the set of # M o yi is a ranking if for all i,j, k, fi, v the following 
hold: 

OiO^oyi y 9^ oyj , 

9» o yi y 9 V o y k 9^ o yi y 9i9 u o y k . 
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If \n\ > \u\ implies # M oyi y- 9 V o y k for all j, k, then the ranking is orderly. If j > k 
implies 9^ o yi y- Q u o y k for all /i, u, then the ranking is elimination. 

Given a ranking >~, a linear polynomial / G TZl \ {0} has the leading term ad o 
iJ 6 0, 8 6 K, where o y J is maximal w.r.t. >- among all 9^ o y k which appear with 
nonzero coefficient in /. lc(/) := a G /C\{0} is the leading coefficient and lm(/) := doyi 
is the leading monomial. 

A ranking acts in TZl as a monomial order. If F C 7?.^ \ {0}, lm(F) will denote the 
set of the leading monomials and \m.j(F) will denote its subset for the indeterminate yf 
Thus, 

m 

lm(F) = UK(F). 

3 Janet-like Grobner bases 

Given a nonzero linear difference ideal X = Id(G) and a ranking y, the ideal generating 
set G = {g h . . . , g s } C K L is a Grobner basis [7] of X if for all / G X n ft L \ {0}: 

3 3 GG, #G0 : lm(/ ) = o lm(#) . (1) 

It follows that / G X \ {0} zs reducible modulo G: 

f^f':=f-\c(f)9o(g/\c(g)), f'el. 

9 

If /' 7^ 0, then it is again reducible modulo G, and, by repeating the reduction, in finitely 
many steps we obtain 

J G 

Similarly, a nonzero polynomial h G TZl, whose terms are reducible (if any) modulo a set 
F C TZl, can be reduced to an irreducible polynomial h, which is said to be in normal 
form modulo F (denotation: h = NF(h, F)). 

In our algorithmic construction of Grobner bases we shall use a restricted set of re- 
ductions called Janet-like (cf. [13]) and defined as follows. 

For a finite set F C TZl \ {0} and a ranking >-, we partition every lnu^F) into subsets 
labeled by do, . . . , di G Z> , (0 < i < n). Here [0]^ := lm^F) and for i > the subset 
[do, • • • , d>i]k is defined as 

{u G lm k (F) | d = 0, dj = deg^u), 1 < j < i}. 
Denote by hi(u,lm.k(F)) the nonnegative integer 

max{deg i (w) | u, v G [d , • • • , - deg^w). 

If hi(u,\m k (F)) > 0, then 0?* such that 

Si := min{degi(u) - deg^w) | u, u e [d , • • • , di-i]*, deg t (w) > deg^w)} 
is called a difference power for / G F with lm(/) = w. 
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Let DP(f, F) be the set of difference powers for / e F, and J~(f, F) := \ 6 be the 
subset of 9 with 

:= {0" | 30" G F) : // - 1/ G Z£ }. 

A Gr6bner basis G of I = Id(G) is called Janet-like [X3J if for all / e I D ^ L \ {0}: 

39GG,^ J^, G) : lm(/) = o lm(^) . (2) 

This implies ^—reductions and the J— normal form NFj(f,F). It is clear that condi- 
tion (j2J) implies condition (fl|). Note, however, that the converse is generally not true. 
Therefore, not every Grobner basis is Janet-like. 

The properties of a Janet-like basis are very similar to those of a Janet basis [14J , but 
the former is generally more compact than the latter. More precisely, let GB be a reduced 
Grobner basis [3], JB be a minimal Janet basis, and J LB be a minimal Janet-like basis 
of the same ideal for the same ranking. Then we have 

Card(GB) < Card(J LB) < Card(J5), (3) 

where Card abbreviates cardinality, that is, the number of elements. 

Whereas the algorithmic characterization of a Grobner basis is zero redundancy of all 
its 5-polynomials O H], the algorithmic characterization of a Janet-like basis G is the 
following condition (cf. [13J ) : 

MgeG^e DP(g, G) : NFj(# o g, G) = . (4) 

This condition is at the root of the algorithmic construction of Janet-like bases as described 
in the next section. 

4 Algorithm for Linear Difference Polynomials 

The following algorithm is an adaptation of the polynomial version [T3] to linear difference 
ideals. It outputs a minimal Janet-like Grobner basis which (if monic, that is, normalized 
by division of each polynomial by its leading coefficient) is uniquely defined by the input 
set F and ranking y. Correctness and termination of the algorithm follow from the proof 
given in [13]; in so doing the displacement of some elements of the intermediate sets G 
into Q at step 13 provides minimality of the output basis. The algorithm terminates when 
the set Q becomes empty in accordance with (j3J). 

The subalgorithm Normal Form(p, G, y) performs the Janet-like reductions (Sec- 
tion [3]) of the input difference polynomial p modulo the set G and outputs the Janet-like 
normal form of p. As long as the intermediate difference polynomial h has a term Janet- 
like reducible modulo G, the elementary reduction of this term is done at step 4. As 
usually in the Grobner bases techniques [3], the reduction terminates after finitely many 
steps due to the properties of the ranking (Section 



4 



Algorithm: Janet-like Grobner Basis(F, y) 



Input: F C TZ L \ {0}, a finite set; y, a ranking 


Output: G, a Janet-like basis of Id(F) 


1 


choose / G F with the lowest lm(/) w.r.t. y 


2 


G:={f} 


3 


Q:=F\G 


4 


while Q ^ do 


5 


h := 


6 


while Q ^ and /i = 0do 


7 


choose p £ Q with the lowest lm(p) w.r.t. >- 


8 


Q-=Q\ M 


9 


ft, := Normal Form(p, G, >-) 


10 


od 


11 


if /i ^ then 


12 


for all g G G such that lm(g) = 8^ o lm(/i), /i > do 


13 


Q:=QU{£}; G:=G\M 


14 


od 


15 


G :=GU{/i} 


16 


Q:=QU{0f>og\geG, 8 13 G DP(g, G) } 


17 


fi 


18 


od 


19 


return G 


Algorithm: Normal Form(p, G, y) 



Input: p G 72. x, \ {0}, a polynomial; G C IZl \ {0}, a finite set; 
a ranking 

Output: h = NFj(p, G), the J — normal form of p modulo G 
1: h := p 

2: while h and ft has a monomial it with nonzero coefficient 

b G K such that w is J— reducible modulo G do 
3: take g G G such that u = 9^ o lm(g) with 6» 7 G J^, G) 
4: ft:= h/b-6"o(g/lc(g)) 
5: od 

6: return ft 



An improved version of the above algorithm can easily be derived from the one for 
the involutive algorithm [H] if one replaces the input involutive division by a Janet-like 
monomial division [13] and then translates the algorithm into linear difference algebra. 
In particular, the improved version includes Buchberger's criteria adjusted to Janet-like 
division and avoids the repeated prolongations 9@ o g at step 16 of the algorithm. 
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5 Computational aspects 



The polynomial version of algorithm Janet-like Grobner Basis is implemented in 
its improved form in C++ [13] as a part of the specialized computer algebra system 
GINV [16]. It has disclosed its high computational efficiency for the standard set of 
benchmarks^. If one compares this algorithm with the involutive one [H] specialized to 
Janet division, then all the computational merits of the latter algorithm are retained, 
namely: 

• Automatic avoidance of some useless reductions. 

• Weakened role of the criteria: even without applying any criteria the algorithm is 
reasonably fast. By contrast, Buchberger's algorithm without applying the criteria 
becomes unpractical even for rather small problems. 

• Smooth growth of intermediate coefficients. 

• Fast search of a polynomial reductor which provides an elementary Janet-like reduc- 
tion of the given term. It should be noted that as well as in the involutive algorithm 
such a reductor, if it exists, is unique. The fast search is based on the special data 
structures called Janet trees U2 



• Natural and effective parallelism. 

Though one needs intensive benchmarking for linear difference systems, we have solid 
grounds to believe that the above listed computational merits hold also for the difference 
case. 

As this takes place, computation of a Janet-like basis is more efficient than computa- 
tion of a Janet basis by the involutive algorithm [TJ]. The inequality for monic bases 
is a consequence of the inclusion 



GB C J LB C JB . (5) 

There are many systems for which the cardinality of a Janet-like basis is much closer to 
that of the reduced Grobner basis than the cardinality of a Janet basis. Certain binomial 
ideals called toric form an important class of such problems. Toric ideals arise in a number 
of problems of algebraic geometry and closely related to integer programming. For this 
class of ideals the cardinality of Janet bases is typically much larger than that of reduced 
Grobner bases [T3]. For illustrative purposes consider a difference analogue of the simple 
toric ideal [131 HH] generated in the ring of difference operators by the following set: 

{ Ox ~ Qy®zi Qx^w ~ 0y, Q x Qy ~ Qz®w } • 

The reduced Grobner basis for the degree-reverse-lexicographic ranking with 9 X y 6 y y 
9 Z y 6 W is given by 



1 Cf . the web page http : // invo . j inr . ru 



6 



The Janet-like basis computed by the above algorithm contains one more element 9l9 w -9\ 
whereas the Janet basis adds another six elements to the Janet-like basis [T3j . 

The presence of extra elements in a Janet basis in comparison with a Janet-like basis is 
obtained because of certain additional algebraic operations. That is why the computation 
of a Janet-like basis is more efficient than the computation of a Janet basis. Both bases, 
however, contain the reduced Grobner basis as the internally fixed [H] subset of the 
output basis^]. Hence, having any of the bases computed, the reduced Grobner basis is 
easily extracted without any extra computational costs. 



6 Nonlinear Difference Polynomials 

In this section we follow the paper |llj and define difference standard bases which gener- 
alize the concept of Grobner bases to arbitrary ideals in the ring TZ = /C-fj/ 1 , . . . , y m } of 
difference polynomials. 

A total ordering y on the set Ai of difference monomials 

M:={(e 1 oy 1 yi...(e rn oy m y™\e j eG, ^gz> , i<j<m} 

is admissible if it extends a ranking and satisfies 
(\/teM\{l}) [tyl] A (V#g0) (Vt,v,weM) [v yw<=^t-0o V y t- 0o W ]. 

As an example of admissible monomial ordering we indicate the lexicographical mono- 
mial ordering compatible with a ranking. 

Given an admissible ordering every nonzero difference polynomial / has the leading 
monomial lm(/) G At with the leading coefficient lc(/). In what follows, every nonzero 
difference polynomial is to be normalized (i.e., monic) by division of the polynomial by 
its leading coefficient. 

If for v,w G Ai the equality w = t ■ 9 o v holds with 9 G O and t G Ai we shall say 
that v divides w and write v \ w. It is easy to see that this divisibility relation yields a 
partial order. 

Given a difference ideal I and an admissible monomial ordering >~, a subset G C X is 
its (difference) standard basis if Id(G) = I and 

(V/eI)(3 5 GG) [lm(s)|lm(/)]. (6) 

As in differential algebra [12], if a standard basis is finite it is called Grobner basis. 

A polynomial p G TZ \ {0} is said to be head reducible modulo q G TZ \ {0} to r if 
r = p — m ■ 9 o q and m G Ai, 9 G are such that \m(p) = m ■ 9 o lm(g). In this case the 
transformation from p to r is an elementary reduction and denoted by p — > r. Given a set 

F C TZ \ {0}, p is head reducible modulo F (denotation: p — >) if there is / G F such that 

F 

p is head reducible modulo /. A polynomial p is head reducible to r modulo F if there is 
a chain of elementary reductions 

p — > pi — > p 2 — > ■ ■ ■ — >■ r . (7) 

F F F F 

2 In the improved versions of the algorithms. 
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Similarly, one can define tail reduction. If r in (EJ) and each of its monomials is neither 
head nor tail reducible modulo F, then we shall say that r is in normal form modulo F 
and write r = NF(p, F). A polynomial set F with more then one element is interreduced 
if 

(V/ GF) [/ = NF(/,F\ {/})]. (8) 

Admissibility of >~, as in commutative algebra, provides termination of chain (J7|) for 
any p and F. In doing so, NF(p, F) can be computed by the difference version of a 
multivariate polynomial division algorithm J2U1 EJj- If G is a standard basis of Id(G), 
then from the above definitions it follows 

/eId(G)<*=>NF(/,G) = 0. 

Thus, if an ideal has a finite standard (Grobner) basis, then its construction solves the 
ideal membership problem as well as in commutative [20] [21] and differential [HI 122] 
algebra. The algorithmic characterization of standard bases, and their construction in 
difference polynomial rings is done in terms of difference ^-polynomials. 

Given an admissible ordering, and monic difference polynomials p and q, the polyno- 
mial 

S(p, q) := mi • 9i o p — m 2 ■ 9 2 ° q 

is called S -polynomial associated to p and q (for p = q we shall say that the S -polynomial 
is associated with p) if 

mi • Q\ o lm(p) = m 2 • #2 o lm(g) 

with coprime m\ ■ 61 and m 2 • #2- 

Algorihmic characterization of standard bases: Given a difference ideal X C 1Z and 
an admissible ordering >~, a set of polynomials G C X is a standard basis of X if and 
only if NF(S(p,q),G) = for all ^-polynomials, associated with polynomials in G. This 
result follows from the above definitions in line with the standard proof of the analogous 
characterization of Grobner bases in commutative algebra [2"0j |2"T] and with the proof of 
similar characterization of for standard bases in differential algebra [19]. 

Let X = Id(F) be a difference ideal generated by a finite set F C 1Z of difference 
polynomials. Then for a fixed admissible monomial ordering the following algorithm 
StandardBasis, if it terminates, returns a standard basis G of X. The subalgorithm 
Interreduce invoked in step 11 performs mutual interreduction of the elements in H 
and returns a set satisfying OH]). 

Algorithm StandardBasis is a difference analogue of the simplest version of Buch- 
berger's algorithm (cf. [191 EH EI]). Its correctness is provided by the above formulated 
algorithmic characterization of standard bases. The algorithm always terminates when 
the input polynomials are linear. If this is not the case, the algorithm may not terminate. 
This means that the do while-loop (steps 2-10) may be infinite as in the differential 
case [IS] [22]. One can improve the algorithm by taking into account Buchberger's criteria 
to avoid some useless zero reductions in step 5. The difference criteria are similar to the 
differential ones [19]. 
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Algorithm: StandardBasis (F, >-) 



Input: F C TZ \ {0}, a finite set of nonzero polynomials; 




y, an admissible monomial ordering 


Output: G, an interreduced standard basis of ld(F) 


1 


G := F 


2 


do 


3 


# := G 


4 


for all S'-polynomials s associated with elements in H do 


5 




6 


if g 7^ then 


7 


G:=GU{ 9 } 


8 


fi 


9 


od 


10 


od while G^H 


11 


G :=Int err educe (G) 


12 


return G 



7 The Maple Package LDA 

The package LDA (abbreviation for Linear Difference Algebrajf] implements the involu- 
tive basis algorithm [TJ] for linear systems of difference equations using Janet division. 
In addition, the package implements a modification of the algorithm oriented towards 
Janet-like division [13] and, thus, computes Janet-like Grobner bases of linear difference 
ideals. 

Tabled] collects the most important commands of LDA. Its main procedure JanetBasis 
converts a given set of difference polynomials into its Janet basis or Janet-like Grobner 
basis form. More precisely, let TZ be the difference ring (cf. Section |2]) of polynomials in the 
variables 9^oy k ) fj, g Z> , k = 1, . . . , m, with coefficients in a difference field K, containing 
Q for which the field operations can be carried out constructively in Maple. We denote 
again by TZl the set of linear polynomials in TZ. Given a finite generating set F C TZl for 
a linear difference ideal I in TZ, JanetBasis computes the minimal Janet (-like Grobner) 
basis J of X w.r.t. a certain monomial order (ranking). The input for JanetBasis consists 
of the left hand sides of a linear system of difference equations in the dependent variables 
y 1 , . . . , y m , e.g., functions of xi, . . . , x n . The difference ring 1Z is specified by the lists 
of independent variables Xi, . . . , x n and dependent variables given to JanetBasis. The 
output is a list containing the Janet (-like Grobner) basis J and the lists of independent 
and dependent variables. 

After J is computed, the involutive/,7— normal form of any element of TZl modulo J 
can be computed using InvReduce. Given p G TZl representing a residue class p of the 
difference residue class ring TZ/X, InvReduce returns the unique representative q G TZl 
of p which is not involutively/,7— reducible modulo J. A /C-basis of the vector space 



j The package LDA is downloadable from the web page http://wwwb.math.rwtri-aachen.de/Janet 
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1Zl/{I H TZl) is returned by ResidueClassBasis as a list if it is finite or is enumerated 
by a formal power series [25] in case it is infinite. For examples of how to apply these two 
commands, cf. Section O 

Given an afflne (i.e. inhomogeneous) linear system of difference equations, a call of 
CompCond after the application of JanetBasis returns a generating set of compatibility 
conditions for the afflne part of the system, i.e. necessary conditions for the right hand 
sides of the inhomogeneous system for solvability. 

Moreover, combinatorial devices to compute the Hilbert series and polynomial and 
function etc. [UJ are included in LDA. 

For the application of LDA to the reduction of Feynman integrals, a couple of spe- 
cial commands were implemented to impose further relations on the master integrals: 
By means of AddRelation an infinite sequence of master integrals parametrized by in- 
determinates which are not contained in the list of independent variables is set to zero. 
Subsequent calls of InvReduce and ResidueClassBasis take these additional relations 
into account (cf. Section E}. 

LDA provides several tools for dealing with difference operators. Difference operators 
represented by polynomials can be applied to (lists of) expressions containing y 1 , . . . , 
y m as functions of the independent variables. Conversely, the difference operators can be 
extracted from systems of difference equations. Leading terms of difference equations can 
be selected. 



Table 1: Main commands of LDA 



JanetBasis 


Compute Janet (-like Grobner) basis 


InvReduce 


Involutive / J — reduction modulo Janet (-like Grobner) basis 


CompCond 


Return compatibility conditions for inhomogeneous system 


HilbertSeries etc. 


Combinatorial devices 


Pol2Shift, Shift2Pol 


Conversion between shift operators and equations 


Some interpretations of commands for the reduction of Feynman integrals: 


ResidueClassBasis 


Enumeration of the master integrals 


AddRelation 


Definition of additional relations for master integrals 


ResidueClassRelations 


Return the relations defined for master integrals 



We consider difference rings containing shift operators which act in one direction only. 
If a linear system of difference equations is given containing functions shifted in both 
directions, then the system needs to be shifted by the maximal negative shift in order 
to obtain a difference system with shifts in one direction only. However, LDA allows to 
change the shift direction globally. 

Unnecessary computations of involutive reductions to zero are avoided using the four 
involutive criteria described in [HJ [231 El] . Fine-tuning is possible by selecting the criteria 
individually. 

The implemented monomial orders/rankings are the (block) degree-reverse-lexicogra- 
phical and the lexicographical one. In the case of more than one dependent variable, 
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priority of comparison can be either given to the difference operators ("term over posi- 
tion" ) or to the dependent variables ( "position over term" / elimination ranking) . 

The ranking is controlled via options given to each command separately. The other op- 
tions described above can be set for the entire LDA session using the command LDAOptions 
which also allows to select Janet or Janet-like division. 

8 Generation of finite difference schemes for PDEs 

We consider the Laplace equation u xx + u yy = and rewrite it as the conservation law 



and using the midpoint integration method we obtain the following discrete system: 



where 9 X and 9 y represent the right-shift operators w.r.t. x and y, e.g., (6 X o u y )(x,y) = 
u y (x + l,y). 

We show how to use LDA to find a finite difference scheme for the Laplace equation: 

> with (LDA) : 

We enter the independent and the dependent variables for the problem (ux > uy > u): 

> ivar: = [x,y]: dvar : = [ux , uy , u] : 

Next, we translate ([9]) into the input format of JanetBasis. Note that one can in general 
use AppShif tOp to apply a difference operator given as a polynomial similar to the ones 
in to a difference polynomial. 

> L:=[2*h*ux(x+l,y)-u(x+2,y)+u(x,y) , 2*h*uy(x,y+l)-u(x,y+2)+u(x,y) , 

> 2*h*(ux(x+2,y+l)-ux(x,y+l))+2*h*(uy(x+l,y+2)-uy(x+l,y))] : 

Then we compute the minimal Janet basis of the linear difference ideal generated by 
L w.r.t. a ranking which compares the dependent variables prior to the corresponding 
difference monomials ("position over term" order; this ranking is chosen when using the 
option 2 as below). The least element of this Janet basis is by construction a difference 
polynomial which does not contain any monomial in ux and uy because ux > uy > u. 

> JanetBasis (L,ivar, dvar, 2) [1] [1] ; 




Adding the integral relations 




-% - e x e 2 y )o Uy + (6% - 6 y )o Ux = 0, 

2Ah6 x o Ux - (91 -l)o U = 0, 
2Ah9yO Uy - {01 - l)ou = 0, 



(9) 
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u(x + 4, y + 2) - 4 u(x + 2, y + 2) + u(x, y + 2) + u(x + 2, y + 4) + u(x + 2, y) 

The computation takes less than one second of time on a Pentium III (1 GHz). 
Dividing this difference polynomial by 4h 2 we obtain the following finite difference 
scheme: 

D*(u jk )+Dl(u jk ) = 0, 

where 

^■(«i*) = 

and 

n 2/„, x U jk+2 ~ 2u jk + U jk „ 2 
D k(Ujk) = ^ 

are discrete approximations of the second order partial derivatives occurring in Laplace's 
equation. 

9 Reduction of Feynman integrals 

In order to demonstrate how to use LDA for the reduction of Feynman integrals, we 
consider a simple one-loop propagator type scalar integral with one massive and another 
massless particle: 

1 f d d s 

f(k,n) := J fcjn 



om-d/2 I pk r>n 
bn J r s—q,m r s,0 

(Here k, n are the exponents of the propagators.) 

The basis integrals for this example and the corresponding reduction formulae were 
found and studied by several authors (see, e.g., [2(i. 27]). Here we apply the Grobner 
basis method, as implemented in LDA, directly to the recurrence relations which have the 
form: 

[d-2k-n- 2m 2 kl+- 
n2+(l" - q 2 + m 2 )} f(k + 1, n + 1) = 0, 

[n- k- kl + {q 2 + m 2 -2-)- { > 

n2+(l- -q 2 + m 2 )} f(k + 1, n + 1) = 0, 

where 

^/(M) = /(^ ± l,n), 2 ± /(A;,n) = /(fc,n± 1). 
In addition, it is known that 

f(k + i,n + j) = Vz<0 Vj (11) 

which we will take into account later. 

> i var : = [k , n] : dvar: = [f]: 

We enter the recurrence relations ( ITU]) : 

> L : =[(d-2*k-n) *f (k+1 ,n+l) -2* m "2*k*f (k+2 ,n+l) -n*f (k,n+2) - 

> n* (m"2-q"2) *f (k+1 , n+2) , (n-k) *f (k+1 , n+1) -k* (q~2+m~2) *f (k+2 , n+1) + 

> k*f (k+2,n)-n*f (k,n+2)-n*(m"2-q~2)*f (k+1, n+2)] : 
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> JanetBasis (L , ivar , dvar) : 

Again, the computation time is less than one second. Now, the master integrals are given 
by: 

> ResidueClassBasis (ivar , dvar) ; 

[f(Jfe, n), f(Jfe, n + 1), f(Jfe + 1, n), f(Jfe, n + 2), f(Jfe + 1, n + 1), f(Jfe + 2, n)] 

(fTTj) implies additional relations on the master integrals. (Here, j is recognized as not 
being contained in ivar and thus serves as a parameter to define the additional relations.) 

> AddRelat ion (f (k,n+j )=0, ivar , dvar) : 
The list of master integrals now becomes: 

> ResidueClassBasis (ivar, dvar) ; 

[f{k + 1, n), f(k + 1, n + 1), f(k + 2, n)] 

Next, we recompute the Janet basis for m = 0: 

> m:=0: J : =JanetBasis (L , ivar ,dvar) : 

For the special case where m = 0, we impose the relation f(k + i,n) = for all i: 

> AddRelat ion (f (k+i ,n)=0, ivar , dvar) : 
Now, we are left with one master integral: 

> ResidueClassBasis (ivar , dvar) ; 

[f(k + l,n + 1)] 

We reduce f(k + 2, n + 3) modulo J taking also the additionally imposed relations on the 
master integrals into account. (Here, the option "F" lets InvReduce return the result in 
factorized form.) 

> InvReduce (f (k+2,n+3) , J, "F") ; 

-((2n + A-d + 2k)(2n + 2-d + 2k)(2k + n-d)(n + 3-d + k) 

(n + 2 - d + k) i(k + 1, n + l))/((n + 1) (2 n - d + 4) n q & k (d - 2 k - 2)) 

Using ResidueClassRelations one can display the relations imposed on the master in- 
tegrals: 

> ResidueClassRelations (ivar, dvar , [i,j]); 

[i(k,n + j), f(k + i, n)\ 

The difference operators occurring in the last result can be extracted as polynomials in 

> Shif t2Pol (% , ivar , dvar , [delta [k] , delta [n] ] ) ; 
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10 Conclusion 



The above presented algorithm Janet-like Grobner Basis is implemented, in its im- 
proved form, in the Maple package LDA, and can be applied for generation of finite 
difference approximations to linear systems of PDEs, to the consistency analysis of such 
approximations [TO] , and to reduction of some loop Feynman integrals. 

Alternatively, the Grobner package in Maple in connection with the Ore algebra pack- 
age [B] can be used to get the same results. 

Two of these three applications were illustrated by rather simple examples. The first 
difference system (discrete Laplace equation and integral relations) contains two inde- 
pendent variables (x, y) and three dependent variables (u,u x ,u y ). The second system 
(recurrence relations for one-loop Feynman integral) also contains two independent vari- 
ables/indices (k,n), but the only dependent variable /. The second system, however, 
is computationally slightly harder than the first one because of explicit dependence of 
the recurrence relations on the indices and three parameters (d,m 2 ,q 2 ) involved in the 
dependence on indices. 

Dependence on index variables and parameters is an attribute of recurrence relations 
for Feynman integrals. Similar dependence may occur in the generation of difference 
schemes for PDEs with variable coefficients containing parameters. Theoretically estab- 
lished exponential and superexponential (depends on the ideal and ordering) complexity 
of constructing polynomial Grobner bases implies that construction of difference Grobner 
bases is at least exponentially hard in the number of independent variables (indices). Be- 
sides, in the presence of parameters the volume of computation grows very rapidly as the 
number of parameters increases. 

The reduction of loop Feynman integrals for more than 3 internal lines with masses is 
computationally hard for the current version of the package. One reason for this is that 
the Maple implementation does not support Janet trees since Maple does not provide 
efficient data structures for trees. 

Another reason is that in the improved version of the algorithm there is still some 
freedom in the selection strategy for elements in Q to be reduced modulo G. Though our 
algorithms are much less sensitive to the selection strategy than Buchberger's algorithm, 
the running time still depends substantially on the selection strategy: mainly because of 
dependence of the intermediate coefficients growth on the selection strategy. To find a 
heuristically good selection strategy one needs to do intensive benchmarking with differ- 
ence systems. In turn, this requires an extensive data base of various benchmarks that, 
unlike polynomial benchmarks, up to now is missing for difference systems. 
For the problem of reduction of multiloop Feynman integrals recently some new reduc- 
tion algorithms have been designed (cf. [28] and references therein) that exploit special 
structure of these integrals and by this reason are computationally much more efficient 
then the universal Grobner bases method. 

The comparison of implementations of polynomial involutive algorithms for Janet 
bases in Maple and in C++ [17] shows that the C++ code is of two or three order faster 
than its Maple counterpart. Together with efficient parallelization of the algorithm this 
gives a real hope for its practical applicability to problems of current interest in reduction 
of loop integrals. 
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Thus, for successful application of the Grobner basis technique to multiloop Feynman 
integrals with masses and to multidimensional PDEs with multiparametric variable coef- 
ficients one has not only to improve our Maple code but also to implement the algorithms 
for computing Janet and/or Janet-like difference bases in C++ as a special module of the 



GINV software [16] available on the web page http : / / invo . j inr . ru . 

As to the algorithm StandardBasis, it has not been yet implemented. Another algo- 
rithmic development also aimed at computation of Grobner bases for systems of nonlinear 
difference polynomials is described in recent paper 
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